Effects of orientation and alignment in high-harmonic generation and above threshold 

ionization 
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When molecules interact with intense light sources of femtosecond or shorter duration the rota- 
tional degrees of freedom are frozen during the response to the strong nonperturbative interaction. 
We show how the frozen degrees of freedom affect the measurable signals in high-harmonic gener- 
ation and above threshold ionization. High-harmonic generation exhibits optical coherence in the 
signal from different orientations of the molecule. For ionization, the contributions from differ- 
ent orientations are added incoherently. The effects are demonstrated for realistic alignment and 
orientation schemes. 
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Currently, intense few-cycle pulses of femtosecond du- 
ration are produced in several laser laboratories world- 
wide. It is of fundamental interest to investigate how 
such pulses interact with quantum systems. Gas phase 
atoms and molecules are ideal test systems for these 
studies. As this research field is maturing, interest- 
ing applications are evolving: tomographic reconstruc- 
tion [if, laser-induced electron diffraction [sj and molec- 
ular clocks [11,111; [Hi 13] being examples under current study. 

Molecules are particularly well-suited for such stud- 
ies since they inherently carry the time scales that 
match those of the new laser sources: nuclei move on 
the femtosecond timescale, electrons on the attosecond 
timescale. Hence, the new ultrashort sources are perfect 
for looking inside molecules and for gaining insight into 
nuclear and electronic dynamics. Conversely, if the dy- 
namics of the system under study is well-understood, in- 
formation about the characteristics of the few-cycle pulse 
can be obtained. For instance, in atoms, where an ac- 
curate description of the electrons can be obtained, the 
carrier-envelope phase difference may be extracted in this 
way @, H, 0| . For molecules, the extra nuclear degrees of 
freedom may be used to control the process of interest. 
As an example, the ability to orient molecules with 
respect to an external axis may be used to enhance the 
high-harmonic yield 



11, 12 



i ay 1 



14] 



Solving the dy- 
namics of the molecule in the strong field is, however, 
much more difficult than the atom due to the extra de- 
grees of freedom. In fact, not even the simplest case 
of H^ interacting with strong IR fields has been solved 
in full 6-dimensional time dependent calculations. For- 
tunately, the interesting prospect of ever shorter pulse 
durations introduces a simplification in the description: 
the timescales of rotation and vibration are often much 
longer than the actual applied pulses themselves, and 
therefore some of these degrees of freedom may be treated 
as frozen during the interaction with the field. In the 
present work, we describe how frozen degrees of freedom 
affect the outcome of an experiment in a non-trivial man- 



ner. We show that the influence of the dynamics and the 
formulation of the theory of measurement depend very 
much on the process considered: we obtain completely 
different behavior for above threshold ionization (ATI) 
and high-harmonic generation (HHG) with respect to the 
coherence in the signal from rotational degrees of free- 
dom. The discussion is exemplified using realistic align- 
ment and orientation schemes. 

We consider HHG and ATI in molecules interacting 
with an ultrashort strong laser pulse. In these experi- 
ments, there are many molecules in the laser focus, but 
the phase space density is low and we can use the single- 
particle density operator p(to) for calculations. At time 
t = to prior to any probe or pump pulse, the molecule 
is in a time-independent thermal state at temperature 
T. By definition, p(t ) = exp (— H/ksT)/Z, with parti- 
tion function Z = Tr[exp(-H/k B T)], H the field-free 
molecular Hamiltonian and ks Boltzmann's constant. 
The initial state is resolved on energy eigenstates \a) 
with energy E a . We concentrate on diatomics where, 
prior to the applied pulses, only the electric and vibra- 
tional ground states are populated. Consequently, the 
energy eigenstate is characterized by the angular mo- 
mentum quantum number, J, and its projection on a 
space fixed axis, M, i.e., a. = (J,M). The discussion is 
straightforwardly generalized to more complicated cases 
and the conclusions remain unaffected. There is no decay 
on the timescales we are considering so propagation is de- 
scribed by a unitary operator U(t): U(t)\a) = \^ a (t)); 
p(t) = U(t)p(t )UHt) = £ a Pc,|*a(t)X*«(*)l. with the 
Boltzmann weight P a = exp (— _E a /fcsT) / 'Z . The evolu- 
tion due to U(t) can contain both alignment pulses and 
a subsequent ultrashort probe pulse producing ATI and 
HHG. 

We separate out the relatively slow rotational move- 
ment of the nuclei to obtain f jAf(r e , R, f2, t) « 
ip(r e ,R, t; ft) 4>jm(Q, t), with R the internuclear distance 
and SI = (9, <f>) the spherical polar solid angle composed 
of the usual polar and azimuthal angles. The variables 
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that enter only parametrically in the wave function are 
put after the semicolon. If we consider the response to a 
femtosecond probe pulse, the rotational degrees of free- 
dom can be considered frozen during the probe pulse cen- 
tered at t p , i.e., the full wave function is approximated 
by 

#j M (r e ,R,Sl,t) « ip(r e ,R,t;tt,t p )cl)jM(Sl,t p ). (1) 

In the case of a prealigning or orienting pump pulse be- 
tween to and t p , 0jA/(f2, t p ) is the rotational wave packet 
evolving from Yjm(Q) at time to. If no pump pulse is 
used then <f) JM {Q,tp) = Y JM (ty- 

Treating HHG first, the complex amplitude for the 
emission of harmonics polarized along the unit vector e 
with frequency u>, is obtained from the Fourier transform 
of the dipole acceleration 



A £ (uj) = e ■ / dt i 



(2) 



E tt f«<*«(t)|dl*e,(t)) 



with (d){t) = Tr p(t)d 

the expectation value of the dipole operator d of the 
molecule. The corresponding power density reads [15j, 
1: 



S e (cj) oc |-A e (a/)| 



(3) 



We note that one can observe interferences in the inten- 
sity S(u>) from incoherent members of the ensemble, i.e., 
members belonging to different ex. This effect of inten- 
sity interferences stemming from adding electric fields is 
known as "polarization beats" to distinguish it from co- 
herent quantum beats fl7l ] . 

We insert the wave function |T]) into the expression for 
(d)(t) and obtain 



(d)(t) 



dQG(n,tp)(d} e ' vlb {t;t p ,n) 



(4) 



with the vibronic dipole (d) e ' vlh (t; t p , fi) 
(<if>(t;tp,n)\de\ip{t;t p ,n)) and 

G(n,t p ) = J2Pj\<PjM(n,tp)\ 2 , 

JM 



(5) 



the angular distribution at time t p . The corresponding 
spectrum is obtained from ©-(O 



S(lo, tp) oc 



J dQG(n,t p )A e e vih (io,n,tp) 



(6) 



with A%' vlh (u, f2, t p ) the complex amplitude for gener- 
ation of high-harmonics at frequency u> and polariza- 
tion e from the electronic dipole (d) e ' vlb (t; O) in a 
molecule fixed at Q. If no alignment pulses are used, 
G(Q, tp) is isotropic. This follows from G(fl,t p ) = 
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FIG. 1: (Color online) Harmonic generation from CO. The 
signal is polarized along the linearly polarized driving laser of 
wave length 800 nm and intensity 2 x 10 14 W/cm 2 . (a) HHG 
yield for (i) an isotropic, unaligned ensemble, (ii) perfect ori- 
entation along the driving laser polarization (||) and (iii) per- 
fect orientation perpendicular to the driving laser polarization 
(_L). (b) HHG yield from CO after field- free orientation by 
a half cycle pump pulse. The delay of the probe pulse with 
respect to the orienting pulse is denoted by t p . 



EjPjEmIYjmWI 2 = Ej^(2J+1)/(4tt) which is 
indeed independent of angles. In this case, the spectrum 
arises from the coherent summation of amplitude contri- 
butions from different orientations, A|' vlb (w, f2, t p ), inde- 
pendent on the temperature of the sample. In general, 
the signals always contain optical coherences, except in 
the trivial case G(fl,t p ) — S(Cl — f2'). Apart from be- 
ing physically well-justified and displaying the effects of 
the frozen molecular degrees of freedom, the relevance of 
([B| is that the numerical propagation during the ultra- 
short pulse is immensely more manageable when the di- 
mensionality is reduced. For molecules with vibrational 
frequencies much smaller than the inverse duration of 
the applied laser pulse, the vibrational coordinates can 
also be treated as fixed, and one can separate out the 
vibrational part of the vibronic wave function in ^ as 
well. For CO with frozen rotations, Fig. [1] shows the 
results of simulations of the emitted high-harmonics of 
the same polarization as the linearly polarized driving 
800 nm probe field with peak intensity 2 x 10 14 W/cm 2 . 
The model used to calculate the harmonic signal resem- 
bles the Lewenstein model. Further details are given 
in [ijj]- Panel (a) shows three cases corresponding to 

(i) uniform orientational distribution, (ii) perfect parallel 
orientation, i.e., the permanent molecular dipolcs (each 
pointing from the O to the less electronegative C) are all 
directed along the polarization vector of the probe pulse 
and (iii) perfect perpendicular orientation. The panel 
clearly shows that the HHG signals depend critically on 
the direction of the orientation. We first discuss cases 

(ii) and (iii) and return to the unaligned case in the end 
of this section. There is no dependence of the azimuthal 
angle in any of the cases, and so we only consider the po- 
lar part of the solid angle. The occurrence and absence 
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of even harmonics in the signal in (ii) and (iii), respec- 
tively, may be understood by recalling that emission of 
harmonic of order N comprises N+l dipole transitions 
(N absorptions of a laser photon, and the emission of one 
high- harmonic photon). In case (ii), only the compo- 
nent || of the dipole operator parallel to the internuclear 
axis is present and this component has a AA = selec- 
tion rule, with A the absolute value of the projection of 
the electronic orbital angular momentum on the internu- 
clear axis. Initially the molecule is in its £ ground state 
and because of AA = it stays in the manifold of £ 
states. The S state, from which the final recombination 
step occurs, is hence accessible by the absorption of both 
an even and odd number of photons and, consequently, 
both even or odd harmonics are produced. Turning to 
case (iii), only the component _L of the dipole operator 
perpendicular to the internuclear axis is present and this 
component has a A A = ±1 selection rule. Consequently, 
only odd harmonics are observed since the II state, from 
which the recombination occurs, can only be reached by 
the absorption of an odd number of photons. In the un- 
aligned case (i), the situation is analyzed by considering 
the transition operator On = H^L"j 1 (j| i cos# + Tj sin#) 
corresponding to emission of a harmonic of order N for a 
molecule with the permanent dipole oriented at an angle 
8 with respect to the polarization vector of the probe laser 
field. In the limits of parallel (9 — 0°) and perpendicular 
(6 = 90°) orientation we retrieve the results discussed 
above. In general the operator On contains even and 
odd powers of cosines and sines. In the unaligned case 
G(180° - 6,t p ) = G{9,t p ), and we see from Eq. © that 
only combinations of the cosines and sines yielding an 
even function on [0, 7r] will survive, i.e., the terms con- 
taining an even number of the |L operator. From the 
selection rules it is, however, clear that the total num- 
ber of A changing transition must be an even number 
and thus the total number of dipole transitions is even, 
explaining why only odd harmonics are emitted in the 
unaligned case. 

In order to simulate a more realistic orientational dis- 
tribution, we present in (b) the prediction of the har- 
monic signal from an ensemble of partially oriented CO 
molecules. To obtain orientation, a half cycle pulse 
(HCP) with amplitude 870 kV/cm and a duration of 0.5 
ps (FWHM) is followed 4.14 ps later by a linearly polar- 
ized laser pulse of 0.5 ps duration and a peak intensity 
of 5 x 10 12 W/cm 2 . In order to model the orientation we 
solve the time-dependent Schrodinger equation for the 
rotational degrees of freedom subject to the rotational 
constant B — 57.9 GHz, dipole moment p = 0.112 D, 
and polarizability volume components otu = 1.925 A 3 
and a± = 1.420 A 3 The initial rotational tem- 

perature is 5 K. The probe field which generates high- 
harmonics is as in panel (a). We plot the harmonic sig- 
nal at two different delays with respect to the peak of 



the HCP. The time delay t p = 5 ps is chosen to illus- 
trate the case with an almost symmetric orientational 
distribution: G(180° -61,5 ps) ~ G(6, 5 ps), and in this 
case the even harmonics are suppressed as expected from 
arguments similar to the ones used in the discussion of 
the unaligned case in panel (a). At longer delays the 
molecules have time to orient obtaining a maximum after 
8.1 ps with (cos#) = -0.11. At this time delay the CO 
dipoles tend to be pointing opposite to the polarization 
vector of the HCP, and this asymmetric distribution al- 
lows for even and odd harmonics of comparable strength. 

We now turn to a discussion of ATI. The funda- 
mental quantity is the probability W(k) for measur- 
ing the momentum k of the outgoing electron. The 
associated measurement projection operator is Vk = 
\ipu)(ipu \ <8> Ir which projects on an electron scattering 
state \ipZ ) with asymptotic momentum k and leaves the 
nuclei unaffected (Ir). Accordingly, W(k) = (Vk) = 
Tr[p(t)Vk] where p(t) is the density matrix of the sys- 
tem. We evaluate the trace in the position-eigenstate ba- 
sis and obtain W(k) = j dr e J dR{r ei R\p(t)P k \r ei R) = 

J2 a p « I dR\ J dr^ (r e )V a (r e ,R, t)\ 2 . Introducing 
the wave functions (fTJ), we obtain: 



W(k,t p )= J dSlG(Q,t p ) J dRR 2 \A(k;R,t p )\ 2 



(7) 



with G(Q, t p ) defined in ([5]) and the complex amplitude 
A(k;R,t p ) = J dr e ip^*(r e )i/j(r e ,R,t;il,t p ) describing 
the transition for frozen rotations at time t p . The time t 
may be any time after the ionizing pulse. In contrast to 
the HHG signal ((6|) which is optically coherent in nuclear 
orientations, the signal for ATI is obtained as an inco- 
herent summation of contributions from different molec- 
ular orientations. As was the case for the HHG pro- 
cess, the calculations involved in the evaluation of ([7]) 
are vastly simplified if the rotational and/or the vibra- 
tional degrees of motion are frozen during the femtosec- 
ond pulse. Figure [2] shows angular distributions of the 
outgoing electron for N2 and O2 for isotropic and field- 
free aligned samples. We obtain the angular distributions 
by integrating ([7]) over the magnitude of the momentum. 
In the present calculation, we approximate the transition 
amplitude in Q by the molecular strong-field approxima- 
tion amplitude [20j . The degrees of alignment are deter- 
mined by the procedure described in 2jJ , and we use the 
following parallel (perpendicular) polarizability volume 
of N 2 : 2.38 A 3 (1.45 A 3 ) and of 2 : 2.3 A 3 (1.1 A 3 ). 
The linearly polarized aligning pulse has peak intensity 
4 x 10 13 W/cm 2 and duration 59 fs (FWHM). The initial 
rotational temperature is 11 K. For the relatively low ini- 
tial rotational temperature assumed here, the maximum 
alignment occurs at the quarter revival period (2.1 ps af- 
ter the pump pulse for N 2 and 2.9 ps delay for 2 ). The 
corresponding values of (cos 2 9) are 0.71 (N 2 ) and 0.74 
(0 2 ). We present results both for alignment preferen- 
tially along the polarization axis of the probe laser and 
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FIG. 2: (Color online) Angular differential ionization rates. 
We use an 800 nm probe laser polarized along the z axis 
for N 2 (a)-(c) at 2 x 10 14 W/cm 2 and 2 (d)-(f) at 1.2 x 
10 14 W/cm 2 . (a) N2 isotropic, unaligned ensemble, (b) N2 
maximally aligned along the z axis, (c) N 2 maximally aligned 
along the x axis, (d) O2 isotropic, unaligned ensemble, (e) 
O2 maximally aligned along the z axis, (f) 2 maximally 
aligned along the x axis. The numbers adjacent to the z axis 
indicate the scale with respect to the randomly aligned cases 
(a) and (d). The plots present an experimentally realizable 
demonstration of the effects of alignment on ATI. 



for alignment along an axis perpendicular to the probe 
polarization. For N2, the only effect of making G(Cl,t p ) 
anisotropic is a change in the overall scaling. This is due 
to the "Eg symmetry of the initial orbital and is in ac- 
cordance with the predictions of tunneling theory which 
favors electron-ejection along the polarization direction 
[22|]. In the case of O2, we find a more pronounced ef- 
fect of the orientation. The change in angular pattern 
reflects the symmetry of the initial Tl g orbital which has 
zero amplitude along and perpendicular to the molecu- 
lar axis: this nodal structure forbids the electron escap- 
ing along the vertical polarization axis when perfectly 
aligned, hence the change from (d) to (e) and (f). 

In conclusion, we have developed the theory of how to 
deal with frozen nuclear degrees of freedom in HHG and 
ATI. The frozen coordinates affect the motion in a non- 
trivial way and coherence issues depend on the degree of 
orientation. Physically, the phase coherence in HHG with 
respect to different orientations may be understood by 
the superposition principle for the electric field generated 
by an ensemble of oscillators with different orientations. 
The incoherence in the ATI signal follows from the fact 
that we in principle, by looking at the nuclear motion 
after the pulse, may infer the orientation of the molecule 
at the instant of ionization. The present work form the 
theoretical basis for future work on molecules interacting 
with strong few-cycle pulses. 
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